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We consider the quantum melting of a two-dimensional flux lattice at temperature T = in the 
"superclean limit." In this regime, we find that vortex motion is dominated by the Magnus force. 
A Lindemann criterion predicts melting when n v /n v > (3, where n v and n p are the areal number 
densities of vortex pancakes and Cooper pairs, and /3 m 0.1. A second criterion is derived by using 
Wigner crystal and Laughlin wave functions for the solid and liquid phases respectively, and setting 
the two energies equal. This gives a melting value similar to the Lindemann result. We discuss 
the numerical value of the melting field at T = for thin layers of low-T c superconductor, such as 
a — MoGe, and single layers of high-T c materials. 
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I. INTRODUCTION. 

Vortices in the layered high-T c materials have remark- 
ably strong thermal fluctuations, which have been exten- 
sively studied [Q. At sufficiently low temperatures, vor- 
tex lines are also expected to be subject to quantum fluc- 
tuations. Quantum effects should manifest themselves in 
the zero-point motion of vortex lines. If these are large 
enough, the flux lattice can melt even at temperature 
T=0. Indeed, many experiments suggest that vortex lat- 
tice melting, both in high-T c materials |^-^| and in low- 
Tc films and multilayers ||, is strongly influenced by 
quantum fluctuations. 

Several authors have already considered the possibil- 
ity of quantum melting in the high-T c superconductors. 
Blatter and Ivlev Q have examined the influence of quan- 
tum fluctuations at finite temperatures. They estimated 
the shift in the melting curve using a Lindemann crite- 
rion, assuming overdamped dynamics combined with a 
Matsubara formalism. Chudnovsky [Q has studied a hy- 
pothetical two-dimensional (2D) quantum vortex liquid 
state at temperature T = 0. Onogi and Doniach |sj com- 
puted the T — melting field for a 2D superconductor 
using quantum Monte Carlo (QMC) techniques without 
dissipative quantum tunneling. By taking into account a 
fictitious magnetic field arising from the Magnus force on 
the vortex pancakes || , they also found strong numerical 
evidence for fractional quantum Hall (FQHE) states in 
the vortex liquid. Such FQHE states had been predicted 
by several authors |l0| , pd[| , but principally in the context 
of Josephson junction arrays. 

In this paper, we describe two simple, quasi-analytical 
models for estimating the conditions for quantum melting 
of a 2D vortex lattice at T = 0, in which the fictitious 
magnetic field is explicitly included. The same model 
might also apply to one layer in a 3D stack of uncoupled 
layers of high-T c material. The first estimate is a simple 
Lindemann criterion. The second involves a simple com- 
parison of internal energies in the crystalline and liquid 



phases. 



II. FICTITIOUS MAGNETIC FIELD AND 
LINDEMANN MELTING CRITERION. 

In our model, the vortex pancakes experience two types 
of forces: those due to other pancakes, and the Magnus 
force arising from the density of Cooper pairs. We ne- 
glect dissipative forces from the "viscous" normal elec- 
tron background, as may be acceptable in the so-called 
superclean limit p^ |. Of the two remaining forces, the 
Magnus force usually dominates (see below). The melt- 
ing Lindemann melting criterion proves independent of 
the vortex mass. By contrast, in the opposite limit where 
the intervortex forces dominate, the melting field depends 
sensitively on the vortex mass Ja] . 

The so-called Magnus force || acting on a single two- 
dimensional ("pancake") vortex, in its rest frame, is 



2e 

7 p = qvhv x zn p = — v x zB e ff. 



(1) 



Here q v = ±1 is the effective charge of the pancake vor- 
tex, h is Planck's constant, v is the pancake velocity, n p 
is the effective areal number density of Cooper pairs (dis- 
cussed further below), B e ff = <&on p is the fictitious field, 
$0 = hc/2e is the flux quantum, and the film is assumed 
perpendicular to the z-axis. 

We now wish to show that the intervortex force is typi- 
cally small compared to the Magnus force. For simplicity, 
consider a superconducting film of thickness d and Lon- 
don penetration depth A. The direct interaction potential 
between two pancakes separated by r is 
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n(r) = 2e dK ( — 



where Aj^ = A 2 jd is the transverse penetration depth, 
e = $q/(167t 2 A 2 ) gives the energy scale of the interac- 
tion per unit length, and Kq{x) is the modified Bessel 
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function of zeroth order. To estimate the effects of the 
vortex- vortex interaction, we assume that the vortices are 
ordered into a triangular lattice, and calculate the change 
in potential energy per vortex, AUharmi due to har- 
monic vibrations about this lattice. After some algebra, 
this extra energy is found to take the form AUharm — 
^4 (|uo — ui| 2 ). Here 1 is a lattice vector of the trian- 
gular lattice, ui is the displacement of the I th vortex from 
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equilibrium, and x(7) = 

where I — and the primes denote differentiation. 

We estimate this energy as follows. First, since the 
vortex-vortex interaction is assumed small, we neglect 
(uo ■ ui). Secondly, in the weak-screening regime where 
the nearest neighbor intervortex distance ag << X±, the 
summation may reasonably be replaced by an integral. 
With these approximations, and using several identities 
for derivatives of Bessel functions, we finally obtain 
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(3) 



where n v = 2/(\/3ag) is the areal vortex density. 

Similarly, for a pancake of mass m v moving in a fic- 
titious field B e ff, the zero-point energy per pancake 
AU mag for a pancake in the lowest Landau level is 



AU mag = -huff, 



(4) 



where = 2eB e ff/(m v c). 

To show that the zero point motion is usually dom- 
inated by B e ff. we will demonstrate that Huj c <C 
faui^f , where cu c is the frequency for zero-point motion 
of the harmonic lattice in the absence of R e ff- Now 
uj c = ^k/m v , where k is the effective spring constant 
of the harmonic lattice. It follows from eq. (3) that 
k = 2eodim v . 

To compare oj c and oj^ , we use the London estimate 
for the penetration depth A 2 (T) = (m p c 2 ) / (4Trq 2 n p D ) = 
(m p c 2 d)/ (4irq 2 n p ), where rip D is the pair density per unit 
volume, m p is the pair mass, and q the pair charge. Then 
a little algebra reveals that u c <C lo & J* provided that 
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where m p is the Cooper pair mass . As will be shown 
below, n v /n p w 0.1 at the melting point. Then inequal- 
ity (9) is satisfied so long as m v /m p <C 20. Now in 
BiSraCaaCusOg+a, the mass of a single pancake vortex, 
assuming a thickness d ~ 10A (appropriate for a sin- 
gle layer of high-T c material) has been estimated as one 
electron mass 0. Thus, in this regime, the inequality is 
satisfied and AUharm <C AU ma g as required. Hence, in 
calculating melting behavior for vortices of this mass, we 
apparently need consider only AU rnag . Our results based 
on considering only AU mag do indeed give n v /n p s» 0.1, 
thereby confirming the self-consistency of our approach. 



We now obtain a simple Lindemann melting criterion, 
assuming that the dominant contribution to zero-point 
vortex motion arises from B e fj. Although lo & J' clearly 
depends on m v , the corresponding zero-point displace- 
ment does not. We calculate this displacement assuming 
the symmetric gauge for the fictitious vector potential, 
A e ff = |B e /y x r. Then in the lowest Landau level, one 
finds 
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independent of vortex mass. 

According to the Lindemann criterion, melting occurs 
the zero-point amplitude is a certain fraction, say a^, of 
a . In most conventional materials, ~ 0.1 — 0.2. Since 
ao = (2<i>o/v / 3i?) 1 ^ 2 , the Lindemann criterion becomes 
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0.07, 



(7) 



using the estimate a 2 L as 0.02. Thus, the Lindemann 
picture predicts quantum melting at T — at a vortex 
density of around 7% of the effective density of Cooper 
pairs per layer. 



III. LAUGHLIN LIQUID VERSUS WIGNER 
CRYSTAL. 

Next, we describe an alternative way of estimating the 
melting temperature in a 2D lattice. We treat the pan- 
cake vortices as bosons, moving in the effective magnetic 
field B e f f . To describe the bosons, we use a Wigner crys- 
tal (WC) wave function in the solid phase, and a properly 
symmetrized Laughlin wave function in the liquid state. 
The melting point is determined by the condition that 
the energies Ewe and Ell of the solid and liquid states 
should be equal. A related approach has been used to 
discuss melting of the 2D electron lattice in a magnetic 
field 0EI 

The WC wave function may be written 



wc 



as n^- 1 ) 



(8) 



Here V'( r ) denotes the zero- momentum single-particle 
wave-function of the lowest Landau level, S is the sym- 
metrization operator, and A is a normalization con- 
stant. We wish to calculate the averaged vortex-vortex 
interaction energy in this state, i. e. Ewe /(2eodS) = 

(*IEi 1 Ei^ U Ko(^^)|*)/25, where S is the 

sample surface area. 

We simplify the calculation by several approximations. 
First, since a 2 , S> (|uo| 2 ), the wave function symmetriza- 
tion is quantitatively unimportant for calculating Ewe- 
Indeed, for large argument, the single-particle wave func- 
tion ip( r ) decays exponentially, and the overlap inte- 
gral between ip(r — L) and ip(r — I2) is almost zero, 
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unless li = I2. In view of this degree of localization, 
Ewe can be expanded in powers of the small ratio 
(|uo| 2 )/Aj_, keeping only the first two terms. The result 

is E wc /(2e dS) = K/2)E M o K o(at) + ""^T. 
where AUharm is given by eq. (3). The fluctuations 
(|u | 2 ) appearing in eq. (3) are, as noted previously, the 
sum of two parts: one due to B e f f and the other to the in- 
tervortex potential. Of these, the former is usually much 
larger, as noted above, and has already been evaluated in 
eq. (6). We substitute this value into eq. (3) and hence 
into the expression for Ewe ln the limit do <C Aj_, 
one can evaluate this sum numerically. The result is 
very well fitted numerically by the form Xa^o ( at) ~ 

n v J d 2 rK (j^) -O.SOOln^Ain^)- 1.437. Collecting all 
these results, we finally obtain 

Ewe n v 



2e dS 



d z rK 



r 



0.25 n v ln 



(Ain„)- (9) 



-0.719 n v + ^. 

Zn p 

For the liquid phase, the wave function symmetry mat- 
ters since the pancakes are delocalized. We use as a trial 
wave function an (unnormalized) wave function of the 
Laughlin form pjjfl: 



LL,m 
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Here zj = xj 



iyj is the position coordinate of the ] th 
pancake, and all lengths are expressed in units of the 
"magnetic length" £ = (<& / (2irB e f f)) 1 / 2 . Since the vor- 
tex pancakes are bosons, to must be an even integer. In 
the Laughlin theory of the fractional quantum Hall effect, 
1 /to is the filling fraction of the first Landau level. 

Laughlin's prescription for obtaining the minimizing 
value of to is readily translated to the present problem, 
in which the role of charges and magnetic field are re- 
versed. The generalized prescription is that the minimiz- 
ing to occurs when the number density n p of vortices of 
the fictitious magnetic field equals to times the number 
density n v of fictitious charges, i. e. m = n p /n v . 

We next calculate the internal energy of the Laughlin 
liquid at various values of m. With a change of scale, 
the vortex-vortex interaction energy of the liquid be- 

comes E LL /(2e dS) = K/2tt) Jd 2 xK ( I7 -2=) ff (x), 

where g(x) is the dimensionless density-density correla- 
tion function for the Laughlin liquid (normalized to unity 
at large x), and x is a dimensionless coordinate defined 
by x = v^Jixn v . Since g(x) differs significantly from unity 
mainly in the region x < 1, it is convenient to decompose 
the interaction energy as follows: 
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(11) 
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where 7 ss 0.577... is Euler's constant and we have used 
the small- a; approximation for Ko(x). 

As noted by Laughlin, the correlation function g(r) 
for the Laughlin liquid state is just that of the 2D one- 
component classical plasma (OCP), in which the particles 
interact logarithmically. The last term on the right is, to 
within a factor, just the internal energy of the OCP. We 
can therefore use standard numerical results for the OCP, 
as obtained by Monte Carlo methods by Caillol et al [[l6) . 
Using the analytical fit of these authors to their own nu- 
merical results for the integral J Q xdxhxx{g{x) — 1), we 
find 



J xdx\nx(j)(x) — 1^ — -ln^TrAj^ri!, 
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/ n \ 0.74 1 / si 

-0.3755 + 0.4400(^J - -ln^Ai^J - -ln(47r) + 



Hence, the energy of the Laughlin liquid can be written 
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\i v ln^X^n^J- (13) 



/ n \ 0.74 
0.720 n v +0.4400 n J — ) 

\2n p J 

Finally, the zero-temperature melting transition is 
defined by the equation Ewe — Ell, or « 

/ \ °- 74 
0-440 ( j^- J , or equivalently 



^ w 0.09. 

n p 



(14) 



We see that this result agrees remarkably well with the 
Lindemann criterion. 
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IV. DISCUSSION. 

We now evaluate these predictions for two materi- 
als of interest, using a simplified approximation for n p . 
As noted by Ao and Thouless, n p is not simply the 
areal density of Cooper pairs per unit area, but that of 
superconducting Cooper pairs - that is, those not pinned 
by lattice disorder. Since it is not obvious how to evaluate 
this quantity, we simply use the London equation to esti- 
mate n p at zero field. To get n p (B), we use the Ginzburg- 
Landau approximation \(B, 0) = A(0, 0)/[l - B/B c2 } 1/2 , 
where B C 2 is the T — upper critical field, and A(B, T) 
is the magnetic field and temperature dependent pene- 
tration depth. The melting condition, from either the 
Lindcmann criterion or from equating solid and liquid 
energies, is n v /n p = /?, where (3 « 0.1. Substituting the 
above expressions into this melting condition, we obtain 



Ba 



Bo + B c 



(15) 



where B = /3m p c 2 d$ /[47rA 2 (0, 0)g 2 ]. 

First, we apply this result to an amorphous MoGe film, 
an extensively studied 2D extreme Type-II superconduc- 
tor. An amorphous Moo.43Geo.57 film of thickness 30A 
has A(0,0) rs 80001 and B c2 « 10 4 G @. Taking 
B w H (a good approximation in the extreme Type- 
II limit), and using (3 = 0.1, we find B «7x 10 4 G, 
and therefore B m /B C 2 ~ 0.8 — 0.9. This is consistent 
with the observations of Ephron et al [fL7f , who find a 
superconducting-insulating transition at around 10 kG, 
quite close to the estimated B C 2. The transition in JlTj is 
undoubtedly not uncomplicated quantum melting, since 
it occurs in highly disordered samples. Indeed, it is un- 
doubtedly better described as a continuous phase tran- 
sition from a vortex glass to a Cooper pair glass Jl8| . 
Nonetheless, it is gratifying that our predicted field, es- 
timated for a clean sample, falls rather close to the ob- 
served transition. 

Of at least equal interest is possible quantum melt- 
ing in high-T c superconductors. Since our model is 
strictly 2D, we consider only a single layer of a high- 
Tc material. The result may conceivably be extrapo- 
lated to the most anisotropic Cu02-based high-T c mate- 
rials, such as BiSr2Ca2Cu20s+x- Assuming d = lOA and 
A(0, 0) = 1400A, we obtain B Q w 1.5 x 10 6 G. Estimating 
B c2 = 3 x 10 6 G, we find B m w 10 6 G. Since T c is smaller 
and A(0, 0) is larger in an underdoped sample, however, 
we may expect B m also to decrease in such materials. 

Finally, we comment on the connection between our 
results and the calculations of Q. While these authors 
find FQHE-like commcnsuration effects in the flux liq- 
uid state, their observed melting scales with m,, as if 
there were no influence of B e ff on B m . Our simplified 
analytical calculations suggest that B e ff may dominate 
the melting behavior for sufficiently light pancake masses 
(m„ <C 40m e ). Presumably, this influence of B e ff would 
show up in QMC studies at sufficiently low values of nx„ . 



To conclude, we have calculated the quantum melting 
criterion for a 2D vortex lattice at T = 0, by comparing 
the internal energies of the vortex solid and vortex fluid 
states in a hypothetical superclean limit. We find that, at 
sufficiently low vortex masses, melting behavior seems to 
be dominated by a fictitious magnetic field acting on the 
vortices and produced by the Cooper pair density. The 
calculated melting field is close to the superconducting- 
insulating transition observed in certain thin films of 
amorphous MoGe, and may be within reach of pulsed 
magnetic fields in some underdoped Cu02-based high- 
T r materials. 
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